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We report on the first Fermi Large Area Telescope (LAT) measurements of the so-called "extra- 
galactic" diffuse 7-ray emission (EGB). This component of the diffuse 7-ray emission is generally 
considered to have an isotropic or nearly isotropic distribution on the sky with diverse contributions 
discussed in the literature. The derivation of the EGB is based on detailed modelling of the bright 
foreground diffuse Galactic 7-ray emission (DGE), the detected LAT sources and the solar 7-ray 
emission. We find the spectrum of the EGB is consistent with a power law with differential spectral 
index 7 = 2.41 ±0.05 and intensity, 7(> 100 MeV) = (1.03 ±0.17) x 10" 5 cm" 2 s" 1 sr" 1 , where the 
error is systematics dominated. Our EGB spectrum is featureless, less intense, and softer than that 
derived from EGRET data. 

PACS numbers: 95.30.Cq,95.55.Ka,95.85.Pw,96.50.sb,98.70.Sa 



Introduction: The high-energy diffuse 7-ray emission is 
dominated by 7-rays produced by cosmic rays (CR) in- 
teracting with the Galactic interstellar gas and radiation 
fields, the so-called diffuse Galactic emission (DGE). A 
much fainter component, commonly designated as "ex- 
tragalactic 7-ray background" (EGB), was first detected 
against the bright DGE foreground by the SAS-2 satel- 
lite [l[ and later confirmed by analysis of the EGRET 



data The EGB by definition has an isotropic sky 
distribution and is considered by many to be the super- 
position of contributions from unresolved extragalactic 
sources including active galactic nuclei, starburst galax- 
ies and 7-ray bursts ([3] and references therein) and 
truly-diffuse emission processes. These diffuse processes 
include the possible signatures of large-scale structure 
formation []! . emission produced by the interactions of 
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ultra- high-energy CRs with relic photons |5(, the annihi- 
lation or decay of dark matter, and many other processes 
(e.g., [3} and references therein). However, the diffuse 7- 
ray emission from inverse Compton (IC) scattering by 
an extended Galactic halo of CR electrons could also be 
attributed to such a component if the size of the halo 
is large enough (i.e., - 25 kpc) @|. In addition, 7-ray 
emission from CRs interacting in populations of small 
solar system bodies Q and the all-sky contribution of IC 
scattering of sol ar p hotons with local CRs can provide 
contributions [8l4l0l|. Hence, an extragalactic origin for 
such a component is not clear, even though we will use 
the abbreviation 'EGB' throughout this paper. 

In this paper, we present analysis and first results for 
the EGB derived from the Fermi Large Area Telescope 
(LAT) [ll| data. Our analysis uses data from the initial 
10 months of the science phase of the mission. Essential 
to this study is an event-level data selection with a higher 
level of background rejection than the standard LAT data 
selections, and improvements to the instrument simula- 
tion. These have been made following extensive on-orbit 
studies of the LAT performance and of charged parti- 
cle backgrounds. Together, these improvements over the 
pre-launch modelling and background rejection allow a 
robust derivation of the spectrum of the EGB that is not 
possible with the standard low-background event selec- 
tion. 

Data selection: The LAT is a pair-conversion telescope 
with a precision tracker and segmented calorimeter, each 
consisting of a 4 x 4 array of 16 modules, a segmented 
anti-coincidence detector (ACD) that covers the tracker 
array, and a programmable trigger and data acquisition 
system. Details of the on-board and ground data pro- 
cessing are given in 11]. 

The LAT ground processing makes use of the pre- 
launch background rejection scheme described in (llj . 
The standard low-background event selection resulting 
from this multivariate analysis, termed "diffuse" class, 
has a Monte Carlo predicted background rate of ~ 0.1 Hz 
when integrated over the full instrument acceptance > 
100 MeV. On-orbit investigations of the residual back- 
ground of misclassified particles in the diffuse event se- 
lection indicated a higher level than predicted from pre- 
launch modelling. To reduce the residual particle back- 
ground further, we developed an event selection com- 
prised of the following four criteria in addition to the 
standard diffuse event classification: 1) events are re- 
quired to have a multivariate- analysis assigned 7-ray 
probability that is higher than the standard diffuse selec- 
tion, with the required probability an increasing function 
with energy instead of a constant value as for diffuse class 
events; 2) the distance of extrapolated reconstructed par- 
ticle tracks from the corners of the ACD must be higher 
than a set minimum value to remove particles that enter 
the LAT in a region where the ACD has a lower than 
average efficiency; 3) the average charge deposit in the 
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FIG. 1: Comparison of (a) LAT on-axis effective area and (b) 
orbit-averaged CR background rate integrated over the FOV 
between the enhanced low-background event selection and the 
standard "diffuse" event selection. 



silicon layers of the tracker is required to be small; 4) 
the reconstructed transverse shower size of events in the 
calorimeter is within a size range expected for electro- 
magnetic showers. The first two criteria assist in reducing 
the overall level of CR background. The second two crite- 
ria provide an additional veto against hadronic showers 
and heavy ions that leak through the standard diffuse 
event classification. In addition to these analysis cuts 
the particle background modelling has been updated to 
be closer to the observed on-orbit charged particle rates. 
Furthermore, the instrument simulation now takes into 
account pile-up and accidental coincidence effects in the 
detector subsystems that were not considered in the def- 
inition of the pre-launch instrument response functions 
(IRFs) [13 ■ 

Figure [T^, shows the on-axis effective area (A e ff) for 
our enhanced low-background and standard diffuse selec- 
tions, respectively. The A e g for the enhanced selection 
is reduced for energies > 300 MeV with a peak value 
~ 0.74 m 2 compared to ~ 0.84 m 2 for diffuse class events. 
The A c ff systematic uncertainties for our enhanced low- 
background selection are of the same magnitude as those 
for the diffuse class events, evaluated by comparing the 
efficiencies of analysis cuts for data and simulation of ob- 
servations of Vela: 10% below 100 MeV, decreasing to 
5% at 560 MeV, and increasing to 20% at 10 GeV and 
above. FigureQJ) shows the orbit-averaged residual back- 
ground rate of our enhanced low-background and stan- 
dard diffuse selection, respectively, determined from our 
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FIG. 2: Comparison of expected and measured orbit-averaged 
event rates for two CR dominated data samples. 



improved simulation. With our enhanced event selection, 
the predicted background rejection is improved by a fac- 
tor 1.3-10. 

We estimate the uncertainty of the residual CR back- 
ground in our simulation by comparing two CR back- 
ground dominated LAT data samples to the predictions. 
The first data sample (A) contains all events passing 
the on-board filters, corresponding to a minimal back- 
ground rejection level. The second data sample (B) con- 
tains events that pass the "source" classification in the 
standard analysis 11] but fail to pass the more strin- 
gent diffuse selection. This second sample corresponds 
to a very high level of background rejection but is still 
dominated by charged particles compared to the stan- 
dard diffuse selection, particularly at intermediate and 
high Galactic latitudes. These particles are from the ex- 
treme tails of the CR distributions that are difficult to 
reject. To both samples we apply the selection cuts 3) 
and 4), described above, to remove the heavy ion and 
hadronic shower backgrounds observed on-orbit that are 
not modeled with sufficient accuracy in the simulation. 
In addition, to reduce the 7-ray fraction in both sam- 
ples only events from Galactic latitudes \b\ > 45° are 
used. The remaining 7-ray "contamination" is negligible 
for sample A due to the overwhelming CR rate. For sam- 
ple B, the contamination is less than 10% above 1 GeV, 
but is almost 30% at 200 MeV, even at high Galactic lat- 
itudes. To remove this 7-ray contamination from sample 
B we use the intensity maps from our fits to the data 
(see below), combined with the IRFs corresponding to 
sample B, to determine the expected 7-ray rate. This is 
subtracted from the observed rate of sample B events. 

Figure [2] compares the orbit- averaged event rates mea- 
sured by the LAT and predicted by our simulation for 
datasets A and B. At the minimal background rejection 
level represented by sample A we find agreement within 
± 20%. This shows that the bulk of the remaining CR 
background is well described by the simulation after re- 
moving the particular class of heavy ion and hadronic 



shower events mentioned above. For sample B, the agree- 
ment is within +50%/— 30%, indicating the uncertainty 
in the description of the extreme tails of the CR distri- 
butions. As these tails are responsible for the limiting 
background in the present analysis, we adopt the results 
for sample B as representative of the uncertainty on the 
expected residual CR background. 

Analysis: We use data taken in the nominal "scan- 
ning" mode from the commencement of scientific op- 
erations in mid-August 2008 to mid- June 2009. The 
data were prepared using the LAT Science Tools pack- 
age, which is available from the Fermi Science Support 
Center 1 . Events satisfying our enhanced low-background 
event selection, coming from zenith angles < 100° (to 
greatly reduce the contribution by Earth albedo 7-rays) 
and incidence angles within 65° of the LAT z-axis (the 
LAT field-of-view) were used. This leaves 19 Ms of total 
observation time in the data set. The energy-dependent 
exposure was calculated using the IRFs corresponding to 
our enhanced low-background event selection described 
above. 

The photon counts and exposure were further pro- 
cessed using the GaD GET package, part of a suite of tools 
we have developed to analyse the DGE [l3| . Gamma-ray 
skymaps were generated using a HEALPix 14| isopix- 
elisation scheme at order 6 with 9 independent energy 
bins from 200 MeV to 102 GeV with GaDGET used to 
simultaneously fit a DGE model, solar 7-ray emission, 
and sources (described below) to the resulting skymaps. 
We only consider the Galactic latitude range |6| > 10° 
in this analysis where the DGE is more than an order of 
magnitude weaker than in the Galactic plane. 

The model used for the large-scale DGE is based on 
the GALPROP code 2 . Recent improvements include use 
of the formalism and corresponding code for pion pro- 
duction in pp-interactions by 



151 . 1 1 61 ] , a complete recal- 



culation of the ISRF [13], updated Hi and H 2 gas maps, 
including corrections to the total gas column density de- 
rived from dust reddening maps [18j an improved linc- 
of-sight integration routine, and the addition of informa- 
tion from our ongoing studies of the DGE with the LAT 
[lil Eoj ]. Cosmic-ray intensities and spectra are calcu- 
lated using a diffusive reacceleration CR transport model 
for a nominal halo size of 4 kpc, with a rigidity depen- 
dent diffusion coefficient that is consistent with available 
CR data for the B/C and 10 Be/ 9 Be ratios, respectively. 
We also consider bounding halo sizes 2 kpc and 10 kpc, 
with corresponding self-consistently derived diffusion co- 
efficients, since the size of the CR halo is one of the princi- 
pal uncertainties in the DGE foreground. The injection 
spectra for CR protons and primary electrons are cho- 



1 http:/ /fermi. gsfc.nasa.gov/ssc/ 

2 http://galprop.stanford.edu, model id 77XvMM7A. 
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FIG. 3: LAT measured 7-ray intensity with fit results for 
|6| > 10° including statistical and systematic errors. Fit re- 
sults by component are given in Table [J Note LAT data are 
dominated by systematic uncertainties for the energy range 
shown in the figure. 



sen to reproduce after propagation the locally measured 
spectra, including the recently reported Fermi LAT CR 
electron spectrum |21| . Gamma-ray emissivities are cal- 
culated using the propagated CR spectra and intensities 
folded with the appropriate target distributions included 
in the GALPROP code: Hi, H2, and H11 gas distributions 
for 7r°-decay and bremsstrahlung, and the ISRF for IC 
scattering. Gamma-ray intensity skymaps are obtained 
by direct line-of-sight integration of the calculated 7-ray 
emissivities. 



For the dominant high latitude components, 
bremssstrahlung and 7r°-decay emission from Hi 
and Hll in the local Galaxy (7.5 kpc < R < 9.5 kpc) 
and IC emission, the intensities are fit to the LAT data 
via scale factors. We use the GALPROP skymaps as 
templates with the component normalisations per energy 
bin as fit parameters. The sub-dominant high-latitude 
DGE components, bremsstrahlung and 7r°-decay from 
H 2 , as well as Hi and Hn outside the local region defined 
above, are taken from GALPROP predictions and do 
not vary in the fit. All sources with test statistic above 
200 (i.e., larger than ~ 14cr) found in the internal LAT 
9-month source list are included with the flux per energy 
band per source as a fit parameter. Weaker sources are 
included with fluxes derived from the LAT catalogue 
analysis. In addition templates for the intensity of the 
7-ray emission from CRs interacting in the solar disk 
and radiation field that take into account the 

relative exposure as the Sun transits the celestial sphere 
are included with their normalisations as fit parameters. 



Results: Figure [3] shows the 7-ray intensity measured 
by the LAT and the fit results for the Galactic latitude 
range |6| > 10°. Table U summarises the numerical val- 
ues and uncertainties, including the intensity values for 
the individually fitted DGE components that are not dis- 
tinguished in figure [3] for clarity. The residual intensity 
obtained after fitting the DGE model components, so- 
lar emission, and sources is the sum of CR background 
and EGB. The simulation is used to estimate the CR 
background and uncertainty, as described earlier. The 
CR background is isotropic when averaged over the data 
taking period in this paper and is subtracted to obtain 
the EGB intensity. Additional figures for different lati- 
tude bands and regions of the sky can be found online 
ft 

Our formal uncertainty on the EGB comes from the 
fit using the nominal model. However, the RMS of the 
residual count fraction between LAT data and our model 
for energies above 200 MeV is 8.2%, when averaged over 
regions of 13.4 deg 2 to ensure sufficient statistics. This 
is larger than the 3.3% value expected solely from statis- 
tical fluctuations. We also see correlation of the resid- 
ual count fraction with structures in the DGE model 
skymaps. This suggests a limitation in the accuracy 
of the description of the DGE model. We investigated 
the uncertainty on the EGB flux related to the DGE 
components by varying the relevant parameters in the 
model and re-evaluating the fits for |6| > 10°. At high 
latitudes, the model parameters principally affecting the 
DGE are: the change of the IC emission with different 
halo sizes and the calculation of the IC emission using 
the anisotropic/isotropic formalism [23| (IC + halo in Ta- 
ble [TJ , variations of the CR source distribution and Xco 
gradient (CR propagation model) , and how assumptions 
used to derive Hi column densities from radio data and 
dust reddening measurements affect the distribution of 
Hi in the local region (Hi column density). To quan- 
tify the uncertainty connected to the residual count frac- 
tion, we used the nominal model and examined the vari- 
ation of the derived EGB when different subregions of 
the |6| > 10° sky are fitted (Subregions of \b\ > 10° sky). 
No single component dominates the uncertainties shown 
in the lower half of Table [J We caution that the uncer- 
tainties for the model components cannot be assumed to 
be independent. Hence, there is no simple relationship 
between the combination of individual components and 
the total formal uncertainty. 

The large statistics allow sub-samples of the total data 
set to be used as a cross check. We repeated our analysis 
for events passing our enhanced selection with 1) differ- 
ent onboard trigger rates and 2) conversions in the thin 
or thick sections of the tracker ll| . The first sub-sample 
ensures that we have properly estimated the residual CR 
background, while the second checks that the small frac- 
tion of misreconstructed Earth albedo events that en- 
ter the LAT in the back section do not affect the result. 
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TABLE I: Fit results and uncertainties for the EGB and other components for \b\ > 10°. 

Intensity integrated over energy band (cm -2 s — 1 sr -1 ) 



J—JllLnJ. y 111 vJJL, V 


0.2-0.4 


0.4—0.8 


0.8—1.6 


1.6-3.2 


3.2—6.4 


6.4-12.8 


1 2 8-25 6 


25 6-51 2 


51 2—1 02 4 


Scs,1g factor 


x 10" ^ 


x 10~^ 


x 10" ^ 


x 10" ^ 


x 10 — ^ 


x 10 — ^ 


x 10 — ® 


x 10 — ^ 


x 10" 10 


EGB 


2.4 ± 0.6 


9.3 ± 1.8 


3.5 ± 0.6 


12.7 ± 2.1 


5.0 ± 1.0 


14.3 ± 4.0 


6.3 ± 1.5 


2.6 ± 0.7 


11.1 ± 2.9 


Galactic diffuse (fit) 


4.9 ± 0.4 


25.9 ± 1.8 


12.6 ± 1.3 


50.7 ± 7.2 


17.0 ± 3.0 


50.0 ± 10 


17.1 ± 3.6 


6.1 ± 1.4 


19.1 ± 5.2 


Galactic diffuse (model) 


5.0 


26.0 


11.5 


43.3 


14.7 


47.9 


15.7 


5.2 


17.0 


IC (fit) 


1.5±0.1 


6.8 ±0.5 


3.5 ±0.4 


16.1 ± 2.3 


6.6 ± 1.2 


23.3 ±4.9 
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Foreground modeling related uncertainty in cm 2 s 1 sr 1 

Hi column density +0.1/-0.3 +0.1/-1.7 +0.1/-0.9 +0.1/-3.6 +0.1/-1.1 +0.1/-2.4 +0.1/-0.9 +0.1/-0.2 +0.1/-1.1 

IC + halo size +0.1/-0.2 +0.1/-0.8 +0.1/-0.5 +0.1/-1.8 +0.1/-0.5 +0.1/-0.7 +0.3/-0.3 +0.4/-0.1 +2.9/-0.5 

CR propagation model +0.1/-0.3 +0.1/-1.1 +0.1/-0.6 +0.1/-0.8 +0.1/-0.3 +0.1/-1.2 +1.4/-0.1 +0.4/-0.1 +3.0/-0.1 

Subrcgions of |6| > 10° sky +0.2/-0.3 ±0.8/-1.5 ±0.4/-0.9 +1.9/-2.1 +0.7/-0.5 ±2.5/-1.9 +1.0/-1.5 ±0.5/-0.3 +2.7/-0.9 
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FIG. 4: EGB intensity derived in this work compared with 
EGRET-derived intensities taken from table 1 in [2J and ta- 
ble 3 in [24| . Our derived spectrum is compatible with a 
simple power-law with index 7 = 2.41 ± 0.05 and intensity 
/(> 100 MeV) = (1.03 ±0.17) x 10" 5 cm" 2 s" 1 sr" 1 where 
the uncertainties are systematics dominated. 



The derived EGB spectrum for these sub-samples is com- 
pletely consistent with that derived from the full data set 
using the same analysis procedure. 

Finally, we note that our analysis also indicates a sig- 
nificant detection of the combined solar disk and ex- 
tended solar IC emission. This finding will be explored 
in more detail in a separate study. 

Discussion: Figure U shows the spectrum of the EGB 
above 200 MeV derived in the present analysis, and from 
EGRET data @, 24 1. Our intensity extrapolated to 



100 MeV based on the power-law fit, /(> 100 MeV) = 
(1.03 ± 0.17) x 10" 5 cm" 2 s" 1 sr" 1 , is significantly 
lower than that obtained from EGRET data: /egret (> 
100 MeV) = (1.45 ±0.05) x 10" 5 cm" 2 s" 1 sr" 1 % Fur- 
thermore, our spectrum is compatible with a featureless 



power law with index 7 — 2.41 ± 0.05. This is sig- 
nificantly softer than the EGRET spectrum with index 
7EGRET = 2.13 ± 0.03 2\. To check that the different 
spectra are not due to the instrumental point-source sen- 
sitivities, we adopt F(> 100 MeV) = 10" 7 cm" 2 s" 1 , 
comparable to the average EGRET sensitivity, and at- 
tribute the flux of all detected LAT sources below this 
threshold to the EGB. We obtain an intensity I re s(> 
100 MeV) = (1.19 ± 0.18) x 10" 5 cm" 2 s" 1 sr" 1 and 
a spectrum compatible with a power-law with index 
7 res = 2.37 ± 0.05. Therefore, the discrepancy cannot 
be attributed to a lower threshold for resolving point 
sources. Our EGB intensity is comparable to that ob- 
tained in the EGRET re-analysis by [24{ with an updated 
DGE model, I S mr(> 100 MeV) = "(1.11 ± 0.1) x 10" 5 
cm" 2 s" 1 sr" 1 . However, our EGB spectrum does not 
show the distinctive harder spectrum above > 1 GeV 
and peak at ~ 3 GeV found in the same EGRET reanal- 
ysis. We note that the LAT-measured spectra are softer 
above > 1 GeV than those measured by EGRET also for 
the DGE at intermediate latitudes [20| and for the Vela 
Pulsar 
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THE SPECTRUM OF THE ISOTROPIC DIFFUSE 
GAMMA-RAY EMISSION DERIVED FROM 
FIRST- YEAR FERMI LARGE AREA 
TELESCOPE DATA. SUPPLEMENTARY ONLINE 
MATERIAL. 

Abstract 

Supplementary material concerning the analysis pre- 
sented in " The Spectrum of the Isotropic Diffuse Gamma- 
Ray Emission Derived From First- Year Fermi Large Area 
Telescope Data" is presented here. 



Galactic diffuse model 

Figure [5] displays the individually fitted contributions 
to the galactic diffuse emission arising from inverse 
Compton emission and interaction of CRs with atomic 
hydrogen in the local Galaxy (7.5 kpc < r < 9.5 kpc) 
via bremsstrahlung and pion decay. These contributions 
are omitted in figure 3 of the published article for visual 
clarity, but numerically available in table 1. The inten- 
sity is averaged over galactic latitudes \b\ > 10°. The 
cosmic-ray contamination is subtracted from the EGB 
component shown here. Its contribution can be found in 
figure 3 of the published article. The errors shown in the 
graphs are the quadratic sum of statistical and system- 
atic errors due to uncertainties in the LAT effective area 
and CR background subtraction. 



Comparisons for different sky regions 

Figure 3 in the published article compares the modeled 
with the measured 7-ray intensity averaged over all galac- 
tic latitudes \b\ > 10°. To illustrate the good agreement 
of the used 7-ray emission model, i.e. galactic diffuse 
emission, EGB, point sources and solar emission, over 
the whole fitted region, it is interesting to display this 
comparison for different sub-regions of the sky. 

Figure [5] shows the intensities of different components 
in the model when averaged over independent Galac- 
tic latitude ranges covering low, mid and high galac- 
tic latitudes, 10° < |b| < 20°, 20° < |6| < 60° and 
|6| > 60°. Figure [7] shows the intensities of different 
components in the model when averaged over different 



hemispheres. The hemispheres considered are centered 
at the North Galactic pole (b > 0°), the South Galac- 
tic pole (6 < 0°), the Galactic center (270° < I < 90°) 
and anticenter (90° < I < 270°). Furthermore, Galactic 
latitudes |6| < 10° are excluded from all hemispheres. 

We emphasize that the intensities of the components 
shown in the figures are from the single gamma-ray 
emission model used in the analysis. In particular, the 
EGB component in each figure is identical because it is 



Galactic diffuse total 
isotropic diffuse 
V77?\ atomic hydrogen (local) 

inverse Compton 
— other Galactic diffuse 
l^s^ sources 

data 
model 




10= 

Energy [MeV] 

FIG. 5: LAT measured 7-ray intensity with fit results for 
|6| > 10° including statistical and systematic errors. Inde- 
pendently fitted components of the galactic diffuse emission 
are shown individually. "Other Galactic diffuse" denotes the 
sum of the H2 and non-local Hi and Hn contributions which 
are included using the model intensities (i.e., no fit). Note 
LAT data are dominated by systematic uncertainties for the 
energy range shown in the figure. Ths intensity of the solar 
emission considered in the model is below the range shown in 
the figure. 



isotropic by construction and thus its average intensity 
does not vary across the sky. The accuracy of the model 
with respect to the different sub-regions of the sky shown 
can be judged by comparing the total predicted model 
intensity (black line in figure [6]) to the measured gamma- 
ray intensities per energy band represented by the data 
points. 
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FIG. 6: LAT measured intensity compared to the 7-ray emission model used in the derivation of the EGB averaged over 
different ranges in Galactic latitude. The regions shown are 10° < |&| < 20° (upper left), 20° < |6| < 60° (upper right) and 
|6| > 60° (lower center). "Other Galactic diffuse" denotes the sum of the H2 and non-local Hi and Hn contributions which are 
included using the model intensities (i.e., no fit). Errors include statistical and systematic errors. LAT data are dominated by 
systematic uncertainties for the energy range displayed in the figure. 
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FIG. 7: LAT measured intensity compared to the 7-ray emission model used in the derivation of the EGB averaged over different 
hemispheres on the sky for Galactic latitudes |&| > 10°. The hemispheres shown are centered at the North Galactic pole (upper 
left), the South Galactic pole (upper right), the Galactic center (lower left) and anti-center (lower right). "Other Galactic 
diffuse" denotes the sum of the H2 and non-local Hi and Hn contributions which are included using the model intensities (i.e., 
no fit). Errors include statistical and systematic errors. LAT data are dominated by systematic uncertainties for the energy 
range displayed in the figure. 



